PROCEEDINGS 

°F SCIENCE 



Wave functions of the nucleon and its parity partner 
from lattice QCD 



Nikolaus Warkentin*", Vladimir M. Braun'', Meinulf Gockeler", Thomas 
Kaltenbrunner", Andreas Scliafer''', Gerrit Schiierhiolz Yoshiifumi Nakamura"^, Dirk 
Pleiter' , Roger Horsley"^, James M. Zanotti"^, Paul E. L. Rakow"^, Hinnerk Stiiben^ 

"Institutfiir Theoretische Physik, Universitdt Regensburg, 
93040 Regensburg, Germany 

^Yukawa Institute for Theoretical Physics, Kyoto University, Japan 

'^Deutsches Elektronen- Synchrotron DESY and John von Neumann Institut fiir Computing NIC, 

15738 Zeuthen, Germany 
'^School of Physics and Astronomy, University of Edinburgh, 

Edinburgh EH9 3JZ, UK 

^Theoretical Physics Division, Department of Mathematical Sciences, University of Liverpool, 

Liverpool L69 3BX, UK 
f Konrad-Zuse-Zentrum filr Informationstechnik Berlin, 

14195 Berlin, Germany 



E-mail: nikolaus . warkent in @ physik . uni -regensburg . de 



QCDSF collaboration 

We compute moments of distribution amplitudes using gauge configurations with two flavors 
of clover fermions from QCDSF/DIK and operators which are optimized with respect to their 
behavior under the lattice symmetries. The knowledge of these quantities helps in understanding 
the internal structure of hadrons and in the analysis of (semi-)exclusive processes. We present 
results for the nucleon distribution amplitude which suggest that the asymmetries (the deviations 
from the asymptotic form) are smaller than indicated by sum rule calculations. Using the same 
approach we were also able to calculate the same quantities for the A^*(1535), the parity partner 
of the nucleon. These results show a stronger deviation from the asymptotic form. 

DESY 08-140 
Edinburgh 2008/33 



The XXVI International Symposium on Lattice Field Theory 
July 14-19 2008 
Williamsburg, Virginia, USA 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Wave functions of the nucleon and its parity partner from lattice QCD 



Nikolaus Warkentin 



1. Introduction 

Distribution amplitudes |[T], ^, ^ describe the structure of hadrons in terms of valence quark 
Fock states at small transverse separation and are required in the calculation of hard (semi)exclusive 
processes. A simple picture is obtained at very large values of the momentum transfer. In this limit 
the process factorizes and, for example, the magnetic Sachs form factor of the nucleon Gm{Q^) can 
be expressed as a convolution of the hard scattering kernel h{xi,yi, Q^) and the leading-twist quark 
distribution amplitude in the nucleon (p{xi,Q^) [||], 

Gm(G') f\dx] f\dy] q>\yi,Q^)h{xi,yi,Q^)q>{xi,Q^), (1.1) 

JO JO 

where [dx] = dx\dx2dx^8{l — LLi-^i)' ^iid is the squared momentum transfer in the hard 
scattering process. However, in the kinematic region IGeV^ <Q^ < lOGeV^, which has attracted 
a lot of interest recently due to the JLAB data [^, ^ for Gm, the situation is more complicated. 
Here calculations are possible, e.g., within the light-cone sum rule approach [^, ^. They indicate 
that higher-twist distribution amplitudes become important while higher Fock states do not play a 
significant role. In any case, the distribution amplitudes are needed as input. 

As advocated in the pioneering work f^], lattice QCD is well suited to calculate such non- 
perturbative quantities. Our recent calculation of moments of the nucleon distribution ampli- 
tudes shows that they can be determined reasonably well on the lattice. Furthermore, using the 
same methods we were able to determine distribution amplitudes of the nucleon parity partner 
A^*(1535) with comparable accuracy. We find that the asymmetry of the leading-twist amplitude of 
the nucleon is smaller than in QCD sum rule calculations, in agreement with light-cone sum rules 



[|10p and quark models []11[], which suggest a less asymmetric form. On the other hand, our re- 
sults for A^*(1535) suggest that the asymmetry for the parity partner of the nucleon is considerably 
enhanced. 



2. Basics 

In this section we work in Minkowski space. The leading-twist distribution amplitudes for 
octet baryons and in particular nucleon distribution amplitudes were introduced within the classical 
framework of hard exclusive processes in ^, |3|] The starting point for baryons is the matrix 
element of a trilocal quark operator, which can be written to leading-twist accuracy in terms of 
three invariant functions V , A and T [[l^]: 



Pexplig A^{a)danfaizi 



Pexp /g / Av{T)dT'' ]gp{z2 



abc 



= T {(P • yC)afi[YsN)yfvV{zi-p) + [p ■ mC)apNyfAA{zrp) + {iOf,vp''C)afi[Y>'YsN)yfTT[zrp)] 



Here a,b,c are color indices, \p) denotes a baryon state with momentum p and /, g, h are quark 
fields. We consider these matrix elements for the space time separation of the quarks on the light 
cone with zt = atz (z^ = 0) and Y.i'^i = 1- On the r.h.s. C is the charge conjugation matrix, f[vAj) 
are normalization constants of the leading-twist distribution amplitudes and is the baryon spinor. 
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In momentum space we have 



V{xi) ^ I V{zi ■ p) flexp (/x,(z, • p)) 



(2.1) 



and similarly for A (a;,) and T{xi). The distribution amplitudes V{xi), A(x;) and T{xi) describe the 
quark distribution inside the baryon as functions of the longitudinal momentum fractions x,. They 
also depend on the renormalization scale. 

The moments of distribution amplitudes V''"" = /o^[dx] x[x!^x^ y(xi,X2,X3) are related to ma- 
trix elements of local operators such as 



-abc 



i'D^^...D^'f{0) (C/)„ 



(2.2) 



by Pltw {0\yf'^^{0)\p) = -fvV''""pPp'p'^p"N^{p), with similar relations for A and T, see, e.g., 
[^]. The multiindex Irnfi with / = Ai . . . A/ and similarly for m and n denotes the Lorentz structure 
given by the covariant derivatives 



dp, —igA^ on the r.h.s. of eq. (22). The indices l,m,n 



(without bars) are the total number of derivatives acting on the first, second and third quark, re- 
spectively. The index p corresponds to the uncontracted Lorentz index of the gamma matrices 
in the operators. The leading-twist projection, Pltw, can be achieved, e.g., by symmetrization in 
Lorentz indices and subtraction of traces. 

Since in the nucleon the two quarks / and g are of the same flavor we have additional relations 
for the moments of the distribution amplitudes. 



V 



linn 



V 



mln 



1 hnn 



( mln 



nlmn 



r^mln 



which can be restored from the combination 

1 







Imn 



Imn 



1 Imn 



+ 2T 



lnm\ 



(2.3) 



(2.4) 



by taking into account the isospin symmetry. Hence we have only one independent distribution 
amplitude. In particular the normalization constants are equal, /n = fv = /a = fr, where /n is the 
nucleon wave function normalization constant defined by the choice = 1. For the parity partner 



1 Imn 



of the nucleon A^*(1535) we have exactly the same relations. The combination (p''"" = V''"" — ^ , 
often used in sum rule calculations, can easily be obtained by the relation (p''™ = 20'"'" — 0"™'. 
Due to momentum conservation we have additional relations between lower and higher moments 
of the distribution amplitudes: 0'"'" = + (pK'n+^h + in particular this imphes 



.110^ 



(2.5) 



In the case of next-to-leading twist distribution amplitudes we restrict ourselves to operators 
without derivatives. Then only two additional constants, and appear [|T^, which determine 
the normalization of the twist-four distribution amplitudes. The combinations Ai = Vi — 4V° and 
A2 = 6{Vi —4T^) are also used in the literature. They describe the coupling to the nucleon of two 
independent nucleon interpolating fields used in the QCD sum rule approach. One of the operators, 
^t:, was suggested in [ [l4] ] and the other, in [15]: 



^t(O) 



.,abc 



{Q)Cfu\Q) X (757pt/^(0))„ ^,(0) 



.^abc 
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Their matrix elements are given by 

{0\^,{0)\p) = X^niNN,, {0\^r{0)\p) = hmNNr. (2.6) 

3. Calculation on the Lattice 

It is straightforward to translate the relevant operators to Euclidean space. However, due to the 
discretization of space-time, the mixing pattern of the operators on the lattice is more complicated 
than in the continuum. It is determined by the transformation behavior of the operators under the 
(spinorial) symmetry group of our hypercubic lattice. As operators belonging to inequivalent irre- 
ducible representations cannot mix, we derive our operators from irreducibly transforming multi- 
plets of three-quark operators [ 16 , 1^] in order to reduce the amount of mixing to a minimum. These 



irreducible multiplets constitute also the basis for the renormalization of our operators, which is 
performed nonperturbatively in an RI-MOM-like scheme. In this procedure, the mixing with "total 
derivatives" is automatically taken into account. The numerical results presented in this work were 
obtained using QCDSF/DIK configurations generated with the standard Wilson gauge action and 
two flavors of nonperturbatively improved Wilson fermions (clover fermions). The gauge coupling 
used was j3 = 5.40, which corresponds to a ^ 0.067 fm via a Sommer parameter of ro = 0.467 fm 



[18, T^. The lattice size was 24^ X 48 with pion masses down to 420MeV. 

In the case of the moments considered in this work we can avoid the particularly nasty mixing 
with lower-dimensional operators completely. Note that the operators 

with different multi-indices plfhh but the same Imn are related to the same moments v''"", A'™ and 
r''"", and we make use of this fact not only in order to minimize the mixing problems but also in 
order to reduce the statistical noise by considering suitable linear combinations. 

For the operators without derivatives, i.e., the matrix elements Ai, A2 and /a?, we have per- 
formed a joint fit of all contributing correlators to obtain the values at the simulated quark masses. 
As these are larger than the physical masses a chiral extrapolation to the physical point is required 
in the end. To the best of our knowledge there are no results from chiral perturbation theory to 
guide this extrapolation. Therefore we have adopted a more phenomenological approach aiming 
at Unear (in m\) fits to our data. It turns out that the ratios fN/m% and Xi/mf^ are particularly well 
suited for this purpose. In order to estimate the systematic error due to our linear extrapolation, we 
also consider a chiral extrapolation including a term quadratic in m\ and take the difference as the 
systematic error. The results in the MS scheme at a scale of 2 GeV are given in Table [l[ Note that 
2A1 ~ — A2 for nucleon, a relation that is expected to hold in the nonrelativistic limit due to Fierz 
identities. However, for A^*(1535) we observe a strong deviation from this relation. 



For the higher moments one can proceed in the same way and the constraint ( |2.5| ) is satisfied 
very well. However, the statistical errors in this approach are too large to allow an accurate deter- 
mination of the (particularly interesting) asymmetries. We achieved smaller errors by calculating 
the ratios T?'™ = where 

51 = + + for/ + m + ?i = l, (3.1) 

52 = 2(0<"i+ 0101 + 0^1°) + 0200 + + for/ + m + n = 2. (3.2) 



4 



Wave functions of the nucleon and its parity partner from lattice QCD 



Nikolaus Warkentin 



These ratios are extrapolated linearly to the physical masses. Again, we extrapolate all results also 
quadratically in m\ and take the difference as an estimate for the systematic error of the chiral 
extrapolation. Requiring that the constraint (2^) be satisfied for the renormalized values we can 
finally extract the moments from the ratios. These results are summarized in Table |l|. In Table ^ we 
compare the moments (p''"" as obtained from 0''"" with some other estimates. From these values 
we see that the asymmetry of the A^*(1535) distribution amplitude is more pronounced compared 
to the nucleon and is mostly driven by the first moments. However, in both cases we have the 
approximate symmetry <p''"" <p'""'. 





Asympt. 


LATW 


LATA'*(1535) 


fN 


lO^peV^] 




3.144(61)(29)(54) 


4.417(1 14)(215)(2) 




■10^ [GeV^] 




38.72(70)(43)(106) 


40.88(110)(778)(57) 


h 


103[GeV2] 




76.23{139)(84)(207) 


208.9(47) (384) (42) 




0100 


1 !^ 0.333 


0.3638(11)(68)(3) 


0.4007(14) (48) (12) 




0010 


1 !^ 0.333 


0.3023(10)(42)(5) 


0.2610(18)(13)(16) 




0001 


1 !^ 0.333 


0.3339(9) (26) (2) 


0.3384(11)(35)(4) 




0011 


1 « 0.143 


0.0724(18)(82)(70) 


0.0706(22) (37) (66) 




0101 


1 « 0.143 


0.1136(17)(32)(21) 


0.1213(23)(21)(16) 




0110 


j « 0.143 


0.0937(16)(3)(38) 


0.0943(24)(17)(38) 




0200 


^ ~ 0.095 


0.1629(28) (7) (68) 


0.1825(35)(4)(56) 




0020 


^ ~ 0.095 


0.1289(27)(37)(51) 


0.0962(37) (68) (119) 




^002 


^ ~ 0.095 


0.1488(32)(77)(73) 


0.1429(34)(29)(75) 



Table 1: Comparison of our lattice results (LAT) as obtained from QCDSF/DIK configurations at )3 = 5.40 
for the nucleon (A^) and A^*(1535) at jx^ — 4GeV^ . The first error is statistical, the second (third) error 
represents the uncertainty due to the chiral extrapolation (renormalization). The systematic errors should be 
considered with due caution, see the text for their determination. 





Asympt. 


QCD-SR 


BK 


BLW 


LAT/V 


LAT A?* (1535) 


fN 


lO' [GeV2] 




5.0(5) 


6.64 


5.0(5) 


3.234(63)(86) 


4.544(117)(223) 




■ lO' [GeV2] 




27(9) 




27(9) 


35.57(65)(136) 


37.55(101)(768) 


A2 


lO'lGeV^] 




54(19) 




54(19) 


70.02(128) (268) 


191.9(44)(391) 




<pioo 


1 f» 0.333 


0.560(60) 


0.38 


0.415 


0.3999(37)(139) 


0.4765(33)(155) 




(pOlO 


i f» 0.333 


0.192(12) 


0.31 


0.285 


0.2986(1 1)(52) 


0.2523(20)(32) 




(pOOl 


{ f» 0.333 


0.229(29) 


0.31 


0.300 


0.3015(32)(106) 


0.2712(41)(136) 




^200 


^ f» 0.143 


0.350(70) 


0.18 


0.212 


0.1816(64)(212) 


0.2274(89) (307) 




^020 


i f» 0.143 


0.084(19) 


0.13 


0.123 


0.1281(32)(106) 


0.0915(45)(224) 




^002 


i fs 0.143 


0.109(19) 


0.13 


0.132 


0.1311(113)(382) 


0.1034(160)(584) 




/II 


^ ~ 0.095 


-0.030(30) 


0.08 


0.053 


0.0613(89)(319) 


0.0398(132)(497) 




ipioi 


^ ~ 0.095 


0.102(12) 


0.10 


0.097 


0.1091(41)(152) 


0.1281(56)(131) 




ipiio 


^ 0.095 


0.090(10) 


0.10 


0.093 


0.1092(67)(219) 


0.1210(101)(304) 



Table 2: Comparison of our lattice results (LAT) for the nucleon and A^*(1535) as obtained from 
QCDSF/DIK configurations at j3 = 5.40 using ^O^", ^H", 0^00 ^nd ^^^o jjj ^able |l] to selected sum 
rule results [|^ (QCDSR) and the phenomenological estimates [|l^ (BLW) and [[llj (BK) at the scale 
11^ = IGeV^. 
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Figure 1: Barycentric plot of the models of the distribution amplitudes for nucleon (a) and A^*(1535) (b) at 
~ IGeV using the central values of the lattice results. The lines of constant xi, xi and xt, are parallel to 
the sides of the triangle labelled by X2, xt, and xj, respectively. 



4. Model for Distribution Amplitudes 

Let us now expand the distribution amplitude in terms of orthogonal polynomials P„y chosen 



such that the one-loop mixing matrix is diagonal |21]: 



,1=0 y=0 



".(Mo) 



Taking N = 2 and calculating the coefficients Cnj{lM)) from an independent subset of the moments 
(jUo = 2GeV), we obtain a model function for the distribution amplitude presented in Fig. 1(a). 



While the (totally symmetric) asymptotic amplitude has a maximum for x\=X2=xt,= 

1 /3, inclusion of the first moments (i.e., choosing N = \) moves this maximum in the case of the 
nucleon to x\ ^ 0.46, X2 ~ 0.27, X3 Rr: 0.27 giving the first quark substantially more momentum 
than the others. The second moments then turn this single maximum into the two local maxima 



in Fig. |l(a)| . This also happens in the case of A'^*(1535). However the influence of the second 
moments on the form of the DA is reduced compared to the nucleon due to the stronger asymmetry 
in the first moments. The approximate symmetry in X2 and X3 seen in both cases is due to the 
approximate symmetry (p^'"" « (p'"'" of our results. It is also seen in QCD sum rule calculations 
for the nucleon as well as in several models such as BLW and BK. Since higher-order polynomials 



have been disregarded. Figs. 1(a) and 1(b) should be interpreted with due caution. 
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